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ABSTRACT 

We present differential stellar population models, which allow improved determinations 
of the ages, iron and a-element abundances of old stellar populations from spectral fitting. 
These new models are calibrated at solar abundances using the predictions from classical, 
semi-empirical stellar population models. We then use the predictive power of fully synthetic 
models to compute predictions for different [Fe/H] and [a/Fe]. We show that these new dif- 
ferential models provide remarkably accurate fits to the integrated optical spectra of the bulge 
globular clusters NGC6528 and NGC6553, and that the inferred [Fe/H] and [a/Fe] agree with 
values derived elsewhere from stellar photometry and spectroscopy. The analysis of a small 
sample of SDSS early-type galaxies further confirms that our a-enhanced models provide a 
better fit to the spectra of massive ellipticals than the solar-scaled ones. Our approach opens 
new opportunities for precision measurements of abundance ratios in galaxies. 

Key words: methods: data analysis, globular clusters: general, galaxies: abundances, galax- 
ies: elliptical 



1 INTRODUCTION 

Stellar population models designed to understand the photospheric 
emission from stars in galaxies and star clusters have reached an 



impressive degree of maturity in the last decade (e.g. 


Vazdekis 


1999| |Brazual & Charlot||2003| |Le Borgne et al.||2004| 


Vazquez 


& Leitherer 2005}. 



A remaining problem is that theoretical stellar spectra are not 
mature enough to reproduce real stars at all spectral types and com- 



plete wavelength range in an absolute sense (e.g. Mart ins & Coelho| 
[20071 jBertone et aT] |2008). The most common solution is to re- 
place in the optical wavelength regime the theoretical stellar spectra 
with spectra derived from empirical spectral libraries. We call such 
models "semi-empirical". A draw-back of these models has always 
been that the stellar spectra are restricted to the abundance patterns 
present in the empirical libraries (dominated by solar neighbour- 
hood stars). So, while semi-empirical models provide an excellent 
basis for solar metallicity and [a/Fe]=0 (e.g. Asari et al. 2007 ; Kol- 
|eva et al.|20 08 1, they are biased at non-solar metallicites. 

To overcome these biases, we here use the predictive power 
of the fully theoretical models from Coelho et al. (2007 hereafter 
C07) to compute relative changes in flux as a function of [Fe/H] 
and [a/Fe]. These are then applied to the solar abundance SSPs 
from semi-empirical models, yielding "differential population syn- 
thesis models" that are at the same time highly accurate in an abso- 
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lute sense and whose predictions for non-solar abundance patterns 
are well-defined and unbiased. The new differential models are ac- 
curate enough to reliably determine the (over-)abundance of the 
a-elements from the integrated spectrum of an old stellar popula- 
tion. Predicting and fitting every pixel in the spectrum, as compared 
to predicting the Lick indices "only" (Trager et al. 2000; Thomas 
|et al.|20 04 1, allows one to (1) use the redundancy of the information 
at different wavelengths, for data of lower S/N, (2) analyse spectra 
at higher resolution (3) explore wavelength ranges not covered by 
Lick indices and (4) free oneself of potentially complicated correc- 
tions for instrumental and velocity dispersion. 



2 COMPUTING DIFFERENTIAL MODELS 

Semi-empirical stellar population models provide a reliable zero- 
point for the spectral evolution of stellar populations of solar metal- 
licity (more precisely [Fe/H]=0 and [a/Fe]=0, called "solar" for the 
remainder of this letter). We therefore rely on them for the abso- 
lute calibration (including their stellar evolutionary tracks). For the 
present letter, we used|Bruzual & Charlot| ( |2003l hereafter BC03), 
Le Borgne et al. ( 2003] hereafter PegHR), and an updated version 
of Vazdekis; H1999| hereafter V09) and we refer to these papers for 
a complete description. 

For the pixel-by-pixel corrections (subsequently called "dif- 
ferentials") as a function of [Fe/H] and [a/Fe] we rely on the purely 
theoretical models computed by C07 and we refer to this paper for 
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all details. The models predict the spectra of SSPs in the age range 
3 to 14 Gyr, in the wavelength range 3000 A to 1 .34 yum, and for the 
abundance values [Fe/H] = -0.5, 0.0 and 0.2 and [a/Fe] = 0.0 and 
0.4 at a constant resolution of full width at half-maximum (FWHM) 
= 1 A. We additionally interpolated the models for [Fe/H] = -0.25 
and [cr/Fe] = 0.2 in order to obtain a finer grid. The interpolation is 
done linearly in log(Flux). This is because the absorption caused by 
a single line of an element species is proportional to the exponen- 
tiated opacity, while the opacity itself depends linearly on element 
abundance. We will refer to this grid of 120 SSPs as the C07 mod- 
els. We make the usual assumption that the fully theoretical models 
correctly predict the flux variation as a function of varying [Fe/H] 
and [a/Fe] (e.g. |Trager et al.|2000||Thomas et al.|2004| >. 

The differential recalibration procedure involves the following 
steps: (1) The SSP spectrum by C07 is convolved with a gaussian 
filter and then resampled to achieve the same resolution and binning 
as the semi-empirical SSP spectra used for calibration. (2) At a 
given age, the differential at each pixel as a function of [Fe/H] and 
at fixed [ff/Fe]=0 is defined as the flux ratio of the SSP flux of a 
given [Fe/H] to the spectrum with [Fe/H] = 0.0. (3) Similarly, the 
ratios between the ff-enhanced fluxes and solar-scaled fluxes, at the 
same [Fe/H], give the differential as a function of [o'/Fe]. (4) All 
semi-empirical models were interpolated in log(Flux) to provide 
the same SSP ages as C07. (5) The last step is then to multiply 
the solar SSP fluxes of the semi-empirical models by the correction 
vectors derived from the fully theoretical models. 

As a result we obtain three sets of differential models (one for 
each calibrator) at the same grid points in age, [Fe/H] and [cf/Fe] as 
the C07 models. We will refer to these as BC03+C07, PegHR+C07, 
and V09+C07, respectively. 

For illustration we show one of the derived differentials, 
namely for [Fe/H]=-0.5 and [a/Fe]=0.0 (Z=0.005), in Figure[T] It 
is clear that the main signal of the variation in [Fe/H] is constrained 
to specific wavelength ranges, such as e.g. the MgB feature. For 
comparison with previous work we also show the differentials of 
the semi-empirical models in Figure [T] We interpolated the semi- 
empirical models in log(Flux) to fall on [Fe/H]=-0.fQ 

The differentials derived from the semi-empirical SSPs are 
significantly bluer than their theoretical counterparts. Insight into 
this problem can be gained by overplotting in Figure [T] the differ- 
ential of the C07 model for the same total metallicity (Z=0.005) 
but with [ff/Fe]=0.4 and therefore [Fe/H]=-0.85 (obtained by ex- 
trapolation in log(Flux)). This new differential shows the same ef- 
fect as the semi-empirical models, only even stronger. Indeed, the 
stars with sub-solar [Fe/H] represented in the stellar libraries are 
generally a-enhanced due to the correlation between [Fe/H] and 
[cf/Fe] that is inherent to the solar neighbourhood (see e.g. Figure 
8 from Bensby et al. 2005 1. Semi-empirical models thus are biased 
to [a/Fe]>0 at subsolar metallicities. The new differential models 
are free of this bias and even predict the change of the SSP spectra 
with [a/Fe]. 

We emphasize that the above differential procedure accounts 
for the effects of abundance variations both on stellar evolution and 
on stellar spectra. It is therefore more accurate than methods based 
on the correction of only stellar spectra ( Prugniel et al. 2007] |Cer-| 
|vantes et al.|2007fr . We note that the difference in the stellar evolu- 
tion prescriptions of the C07 and the semi-empirical models should 



1 The models are generally defined in total metallicity Z and not in [Fe/H] ; 
under the assumption of solar abundance ratios, i.e. [a/Fe]=0, [Fe/H]=-0.5 
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Figure 1. The ratios between SSPs of lOGyr at subsolar metal- 
licity (Z=0.005, [Fe/H]=-0.5, [a/Fe]=0) and solar metallicity 
( [Fe/H] =0, [o'/Fe] =0)) for different stellar population models (black 
C07, green BC03, blue PegHR and red V09). The dotted line is C07 at 
the same Z=0.005, but different abundance ratios ([Fe/H]=-0.85 and 
[Q-/Fe]=0.4). All models have been smoothed to a common resolution of 
<x=lA in the upper panel and <x=5A in the lower panel. 



have only a minor influence on the differentials, as illustrated by the 
similarity of the results obtained for semi-empirical models with 3 
different prescriptions in Figure[T] 



3 PIXEL-FITTING OF SPECTRA 

In order to test the applicability of the model spectra on observa- 
tional data, we use a pixel-fitting software previously used and fully 
described in Wa lcher et al.| < [2006] l, hereafter sedfit. In brief, the 
software uses an initial template to obtain a reasonable value for the 
velocity and velocity dispersion of the object through a^ 2 -biased 
random walk. Then, a non-negative least squares algorithm (Law- 
|son & Hanson 1974) performs the inversion onto the stellar popula- 
tion model. Finally, the ^-biased random walk is repeated to yield 
a best estimate for the kinematic parameters. 

Like the many other pixel-fitting softwares in existence (e.g. 
Tojeiro et al.|2007[[Ocvirk et al.|2006]|Cid Fernandes et al.|2005[ 



Koleva et al. 2008 1, this approach has the advantage of being able to 



fit the kinematics and composite stellar populations simultaneously. 
It also provides for a convenient way to mask spectral regions either 
in the observed frame (CCD defects, sky lines) and/or in the rest- 
frame (emission lines, regions badly fit by the models). 

We caution that the fit results are sensitive to the continuum 
normalization method. Continuum-normalized spectral fitting has 
been favoured by e.g. Wolf et al. (2007) over other methods as the 
most reliable to simultaneously retrieve age and metallicity from 
the integrated spectra of stellar clusters. After extensive testing, we 
settled for a simple running mean with exclusion of 3cr outliers as 
the most robust procedure. The wavelength window o ver which the 
mean was evaluated was set to cover 90 A (compare Koleva et al. 



is equivalent to Z=0.005. 



2008). 
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Table 1. Properties of GCs from colour-magnitude diagrams and high- 
resolution stellar spectroscopy 



[tbp] 



Name A 


?e[Gyr] 


[Fe/H] [ar/Fe] 


NGC6528 11 ±2' 
NGC6553 11±2 3 


-0.10 ±0.2 2 +0.1±0.1 2 
-0.20 ±0.1 4 +0.25 ±0.1 3 


Feltzing & Johnson 


j2002|; 2 


Zoccali et al.|j2004|; 3 |Alves-Brito et al.| 




2006"l; Melendez et al. 1 2003 1 



4 APPLICATION TO GLOBULAR CLUSTERS 

We now test the newly prepared differential models on high qual- 
ity globular cluster spectra of NGC6528 and NGC6553 taken from 
Schiavon et al.| ((2005 )0 Only these two bulge GCs have [Fe/H] in 
the range probed by the differential models. For NGC6528 six spec- 
tra are available for consistency checks. Table[T]lists the published 
ages and abundances of these clusters. Not all literature measure- 
ments have assigned errorbars and we have been forced to assign 
hopefully realistic ones ad hoc. 

We fit the spectra of the GCs using the three different differ- 
ential model sets (BC03+C07, V09+C07, PegHR+C07). While an 
extensive test of the C07 models is beyond the scope of this Let- 
ter, we also quote the results obtained with these models alone to 
illustrate the improvement provided by the differential models. The 
best-fitting SSP is determined by simply computing x 2 f° r every 
SSP with the three different parameters age, [Fe/H], [a/Fe]. The 
only purpose of using sedfit in this context is to automatically 
adapt for small shifts in velocity and for the resolution of the data. 

An example fit is shown in Figure [2] Over the limited wave- 
length range of this fit, the model looks remarkable. A more 
quantitative measure is the reduced x 2 , which is between 2 (for 
NGC6528_b_l with PegHR+C07) and 10 (for NGC6528_a_l with 
BC03+C07). While thus the fit is not perfect, we note that an SSP 
with solar abundance ratios yields a reduced^ 2 that is between 10% 
and 50% worse. Combining the realism of semi-empirical models 
and the predictive power of purely theoretical ones into the differ- 
ential models thus leads to a significant improvement in the fit. 

We repeat the fits for different ranges in wavelength. From the 
outset, we exclude the wavelength region below 4000 A, as there 
are still not fully explained problems with semi-empirical models 
in this range (e.g. [Wild et al.|2007)|Walcher et al.|2008) . We also 
exclude the region above 5800 A, which has some dubious features 
in the data. 

The wavelength regions we tested are, in order of increasing 
coverage: 

Rl: 5142-5364 A (MgB and Fe5335 indices) 

R2: 4828-5364 A (main Lick index range) 

R3: 4828-5600 A (red extension) 

R4: 4500-5600 A (blue extension) 

R5: 4000-5800 A (full wavelength range) 

In Figure [3] we show the dependence of the quartet reduced 
X 1 , age, [Fe/H], and [oyFe] on the adopted wavelength range and 
for the spectrum NGC6528_a_2. Based on all plots for the 7 differ- 
ent GC spectra that we used, we find that for wavelength region R2 
(4828 to 5364 A) the parameters [Fe/H] and [ar/Fe] of the best-fit 
SSP are in good agreement with the literature values for all dif- 
ferential models and all observed GC spectra. On the other hand, 
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Figure 2. The spectrum of NGC6528 (black), labelled a_2, as taken from 
| Schiavon et~aT1j2005) and the best-fitting model from PegHR+C07 (red). 
The blue line is the continuum that was used to normalize the observed 
spectrum. The model was fit to the data only in the wavelength range 4828- 
5364 A. 



the results for the other wavelength ranges scatter more and seem 
in particular affected by the age- [Fe/H] degeneracy. In the case of 
Rl, the available signal may not be strong enough. In the cases of 
large wavelength ranges, we expect that all models become less re- 
liable, so small uncertainties in the modelling might dominate the 
fit results. Concerning age, we confirm the results of Koleva et al. 
(2008 ) that the spectral fitting ages for these clusters tend to be 
younger than the CMD ages. This is potentially attributable to con- 
tamination by blue bulge stars or blue straggler stars in the aperture 
covered by the integrated GC spectra and visible in the CMD (G. 
Beccari, priv. comm.) or to dispersion in the Helium abundances 
( |Catelan et al.|2009) , but we note also that the CMD ages do differ 
by a few Gyrs in the literature. 

It seems to be of rather little relevance, which semi-empirical 
model one uses, as inside the regions in wavelength and parame- 
ter space that we test here, the differential models behave similarly. 
We nevertheless note that the PegHR+C07 model yields the lowest 
overall^ 2 but also systematically lower ages (and higher [Fe/H]), 
while fitting with the BC03+C07 models seems to be most stable 
with respect to wavelength range between 4500-5600 A. As ex- 
pected, the fully theoretical models show a larger^ 2 and the derived 
values tend to be less stable than those of the differential models. 

In Figure|4]we show the distribution of^ 2 as a dependence of 
the three parameters age, [Fe/H], and [a/Fe] for the NGC6553 spec- 
trum. The known age- [Fe/H] degeneracy is clearly visible. There is 
also a slight degeneracy between [Fe/H] and [c/Fe], following a 
line of constant total metallicity. 



5 TEST ON SDSS EARLY-TYPE GALAXIES 

We now fit the spectra of a small sample of early-type galaxies 
from the SDSS survey (Abaza jian & Sloan Digital Sky Survey | 
2008 ). Early-type galaxies are known to be enhanced in tf-elements 
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Figure 3. The dependence of the best fit parameters on the wavelength 
range that was used for the fit and on the model. The numbers of the 
wavelength ranges are given in the text of Section [4] We show results 
for one spectrum (NGC6528_a_2) and for all models, i.e. V09+C07 (red), 
PegHR+C07 (bue), BC03+C07 (green), C07 (black). 



(Wort hey et aI1[T992) |Trager et al.]|2000{ |Vazdekis et aT]|200T) 
|Thomas et al.|200*4 1 and we expect this to be reflected in the best 
fit results. We select a subsample of 150 galaxies from the SDSS 
sample of Gallazzi et al. (20061. The objects have been chosen 
in three bins of stellar mass, log(M»/M Q ) = [10-10.3], [10.7-11], 
[1 1.3-1 1.6]. In each mass bin we include only the first 50 galaxies 
with the highest S/N and with A(MgB/<Fe>) within 0.1 of the typ- 
ical value for the mass range according to the linear fit in |Gallazzi| 
|etaL]((2006) (A(MgB/<Fe>) is an empirical estimate of [a/Fe]). 

In order to test the validity of the differential models on the 
SDSS early-type spectra we use sedfit to fit the star formation 
history of the galaxies in a non-parametric way with two sets of 
templates ranging in age from 3 to 12 Gyr: once with the templates 
[Fe/H] = -0.25 and [a/Fe] = and once with [Fe/H] = -0.25 and 
[a/Fe] = 0.2. The fit was carried out using the V09+C07 model in 
the wavelength range 4828 to 5364 A. 

Figure [5] shows that, while low-mass early-types are equally 
well fitted by scaled-solar than by alpha-enhanced models, for more 
massive early-types [a/Fe] =0.2 models provide on average a better 
fit. This reflects the expected mass-dependence of a-enhancement 
in local early-type galaxies. We hold this result to show that the 
differential models represent the spectra of early-type galaxies in 
the SDSS in a realistic way, over the whole spectrum between 4828 
and 5364 A. 

The choice of values for [Fe/H] and [a/Fe] in Figure [5] is mo- 
tivated by the distribution of the sample in these parameters as 
measured using the differential models. Nevertheless, due to the 
degeneracy between [Fe/H] and [a/Fe] noted in Figure 4, simi- 
larly good fits are obtained when increasing [Fe/H] to 0.0, but the 
best fit [a/Fe] would be shifted down by roughly 0.1 in compen- 
sation. A more thorough discussion of the impact of this degener- 
acy, a comparison with the global metallicities derived by Gallazzi 
et al. (2006) and the implications of the results for the [Fe/H] and 
[a/Fe] distribution of SDSS early-type galaxies will be presented 
in a forthcoming paper (Gallazzi et al., in prep.). 



Figure 4. The distribution of x 2 for the NGC6553_a_l spectrum and de- 
termined over the wavelength region 4828-5364 A, depending on the SSP 
parameters age, [Fe/H], and [cr/Fe] using the V09+C07 model. The thicker 
the line the higher the value of^ 2 , in steps of Ax 2 = 1. The squares with 
errorbars show the literature values from Table [T] The upper right panel 
shows that it is important to obtain both [Fe/H] and [a/Fe] correctly, in or- 
der to avoid a biased age measurement. In the lower right panel the locus of 
constant total metallicity Z=0.()17 is indicated by a thin straight line. 



6 SUMMARY 

We have computed new, differential stellar population models. 
These combine (1) the realism of semi-empirical models at solar 
metallicities, due to the recently collected extensive libraries of 
stellar spectra, and (2) the predictive power of purely theoretical 
models, due to recent advances in the computation of purely theo- 
retical stellar spectra ( |Percival et al.|2009[|Coelho et al.|2007fr . 

We fitted the models to high-quality, public spectra of two 
globular clusters in the range of parameters covered by the new 
models. We showed that our best fit age, [Fe/H] and [a/Fe] are 
in excellent agreement with those determined by other methods in 
the literature. Finally we fitted a sample of 150 early-type galaxies 
from the SDSS and show that the reduced x 2 is better for a model 
with [a/Fe] = 0.2 than for one with [a/Fe] = 0.0 for galaxies more 
massive than ~ 10 10 5 MO, again in accordance with independent 
measurements from indices. 

The reduced x 2 values of the models on SDSS spectra are 
close to 1. We caution that this is not yet true for the highest S/N 
data of galactic globular clusters. Nevertheless, the uncertainties in 
the modelling are now smaller than the effects of the varying abun- 
dances over a significant wavelength range (4828 to 5364 A). It 
needs to be pointed out though, that metallicity estimations from 
fitting the full optical spectrum, as done e.g. in jPanter et al.|{2008] > 
and MacArthur et al. (2009), are still rather uncertain. On the other 
hand, spectral fitting allows to use the redundant information on a 
large part of the spectrum and to adapt automatically for changing 
velocity dispersions, thus it yields tight constraints even for a large 
sample of intermediate S/N spectra such as those of the SDSS. 

We conclude that the differential models can reproduce ac- 
curately the detailed spectra of galaxies, thus allowing to measure 
[Fe/H] and [a/Fe] in galaxies on an absolute scale and with high 
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Figure 5. The reduced x 1 of a composite age fit to the spectra of a sam- 
ple of 150 early-type galaxies from the SDSS at fixed [Fe/H]=-0.25 and 
[o-/Fe]=0.2 (abscissa) and [o/Fe]=0.0 (ordinate). Objects are split by mass, 
with the lowest in blue, intermediate in black and high-mass objects in red. 



accuracy. With future improvements in semi-empirical and in the- 
oretical models, differential stellar population models will become 
even more powerful. 
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